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Abstract 



The disorder averaged single-particle Green's function of electrons subject 

to a time-dependent random potential with long-range spatial correlations 

is calculated by means of bosonization in arbitrary dimensions. For static 

disorder our method is equivalent with conventional perturbation theory based 

on the lowest order Born approximation. For dynamic disorder, however, we 

obtain a new non-perturbative expression for the average Green's function. 

Bosonization also provides a solid microscopic basis for the description of 

the quantum dynamics of an interacting many-body system via an effective 

stochastic model with Gaussian probability distribution. 
PACS numbers: 72.10.-d, 05.30Fk, ll.lO.Ef, 71.27.-Fa 
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Over the past 20 years the motion of a quantum particle in a time-dependent random 
potential has been the subject of many works Most authors have focused on the 

case of a single isolated electron, where numerical- as well as non-perturbative analytical 
methods are available . The equally important problem of calculating the average single- 
particle Green's function in the presence of a filled Fermi sea has not received much attention 
0. The purpose of the present paper is to show that the functional integral formulation 
of bosonization offers a new non-perturbative approach to this problem in arbitrary 
dimensions d. 

The bosonization approach in > 1 has recently been re-discovered by Haldane [§, 
who generalized earlier ideas due to Luther [^]. For a detailed description of the geometric 
construction leading to higher-dimensional bosonization we would like to refer the reader 
to Ref. P). Here we briefly recall the basic features of Haldane's construction. The first 
step is the subdivision of the Fermi surface into disjoint patches K"^ and the introduction 
of a collection of local coordinate systems with origins on the Fermi surface (this is called 
an atlas [0). In d dimensions each patch covers an area A'^^^ of the Fermi surface, where 
the cutoff A should be chosen such that within a given patch the curvature of the Fermi 
surface can be locally neglected. For spherical Fermi surfaces, this means that A should 
be chosen small compared with the Fermi wave-vector kp- Each patch is then extended 
into a dimensional box iC^;^ with radial height A and volume A'^^^A. One proceeds by 
defining local density operators associated with the boxes, = Z^k 0"(k)c|^Ck+q, where Ck 
is the annihilation operator of an electron with wave-vector k (for simplicity we consider 
spinless electrons in this work), and the cutoff function B"(k) is unity for k G K"^^ and 
vanishes otherwise. In the functional bosonization approach 0-0 the calculation of the 
interacting Green's function is then mapped via a Hubbard- Stratonowich transformation 
onto the problem of calculating the average Green's function of an effective non-interacting 
system in a dynamic random potential 0", which couples to the local densities p°. Because in 
the high-density limit the probability distribution of the 0"-field is Gaussian, the calculation 
of the average Green's function reduces to the calculation of a trivial Debye- Waller factor. In 
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the field theory hterature 0] the auxihary field 0" is called disorder field. This terminology 
suggests that the functional integral formulation of bosonization can also be used to calculate 
the average Green's function of non-interacting electrons in a dynamic random potential. 
This connection between bosonization and disordered systems is obscured in the operator 
approach |TI|JT^, but it is quite clear in the functional formulation. Indeed, in this work 



we shall show that the functional bosonization approach developed in Refs. provides 
a useful method to calculate the average single-particle Green's function of electrons in a 
stochastic medium. 

We consider non-interacting spinless electrons at high densities subject to an imaginary- 
time random potential f/(r, r). The imaginary-time Green's function Q{t, r', r, r') is defined 
as usual, 
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g{r, r', r, r') 
S(r-v')6*{r-r') , (1) 



where m is the mass of the electrons, fi is the average chemical potential, and 6*{t) = 
^Z^n^"""^"^- Here f3 is the inverse temperature and Un = 27r(n + |)//? are the fermionic 
Matsubara frequencies. We assume that the random potential has a Gaussian probability 



distribution with zero average and general covariance U{r,T)U{r' ,t') = C{r — r',r — r'), 
where the overbar denotes averaging with respect to the probability distribution of the 



random potential. In Fourier space we have then UqU^q = pVCg, where V is the volume of 
the system, and the Fourier components are defined by 

Cg = J^^ drj rfre-'('i-'-"™")C(r, r) , (2) 

and analogously for Ug. Here q = [c[,iujm], and Um = 2nm/i3 are the bosonic Matsubara 
frequencies. 



We are interested in the average Green's function G{r — r', t — t') = Q(r, r', r, r'). Within 
our bosonization approach G{r — r', r — r') is calculated in the most direct way: First we 
obtain the exact Green's function Q{r, r', r, r') for a given realization of the random potential, 



and then this expression is averaged. As explained in detail in Ref. 0], the first step can 
only be carried out if there exists a cutoff <^ kp such that Uq becomes negligibly small for 
|q| ^ Qc- But this means that we have to assume that for wave- vectors |q| ^ qc the Fourier 
coefficients Cq of the covariance function can be neglected. In other words, we can only 
study random potentials with sufficiently long-range spatial correlations, such that qc <^ kp- 
Evidently the most popular model of 5-function correlated disorder cannot be treated within 
our bosonization approach. However, in view of the fact that a random potential with a 
finite correlation range q~^ is expected to lead for distances |r| ^ q^^ to qualitatively 
identical results for single-particle properties as a 5-function correlated random potential, 
this restriction seems not to be very serious. By choosing the above mentioned box cutoffs 
such that gc A, A <^ kp, we obtain cutoff-independent results for physical correlation 
functions at distances |r| ^ q^^. 

Due to the subdivision of the Fermi surface into small patches it is possible to linearize 
the energy dispersion locally within a given patch, — /i]k=k"+q ~ ■ q, where k" is the 
center of patch K'^, and v" is the local Fermi velocity. We then replace Eq.(0) by a first 
order partial differential equation for the patch Green's function ^"(r, r', r, r'). 



-dr + w° ■ Vr - f/(r, r)] g'^iv, r', r, r') = 

5(r-r')5*(r-r') 



(3) 



Because this differential equation is first order and linear, it can be solved exactly by means 



of a trivial generalization of Schwinger's ansatz |]13[. The result can be written as 



^"(r,r',r,r')=G^(r-r',r-r') 



X exp 



V" ■ q 



with 



(4) 
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where 5'^'''~^\r^) is a Dirac (5-function of the d — 1 components of r that are orthogonal 
to v". The Gaussian average of Eq.(|l) yields the usual Debye- Waller factor. 



g^{r, r', r, r') = G"(r - r', r - r') 



with 



The total disorder averaged Matsubara Green's function can then be written as 

X G^{r,T)e^^^si-'^^ . (8) 

Because G'Q(r,r) is proportional to 6'^~^{fJ), we may replace Q^j^(r, r) ^ <5djs('^"'^"' ''") 
Eq.(H), where r" = ■ r, with v° = This completes the solution of the non-interacting 
problem. 

Because disorder and interactions are treated on equal footing in our bosonization ap- 
proach, it is easy to see that interactions simply give rise to another additive contribution 
Q°'{r,T) to the Debye- Waller factor. Hence, for disordered interacting electrons Eq.(^ 
should be modified by replacing Q2is{r, r) by 

Qto*(r,r) = gL(r,r) + gL(r,r) , (9) 

where for a general (possibly retarded) density- density interaction fq the contribution 
Qf„f(r,r) can be written as 

fr r) - — V ./?PA l-cos(q-r-^^r) 

f3V^^' (z^^-v--q)2 ■ ^^^^ 

Here f^^"^ = fq[i + ^o{<l)fq]~^ is the random-phase approximation for the effective inter- 
action, and Ho(g) is the non-interacting polarization. In d = 1 a result similar to Eqs.(^ 
and (|T0|) has also been obtained by Hu and Das Sarma [|T^ by means of a very different 
(and perhaps less elegant) method. Evidently Eq.(^ does not contain interference terms 
describing weak localization effects. These are known to play an important role in the low- 
energy behavior of the average Green's function of an interacting disordered Fermi system in 



the metallic regime. Note that diagrammatically weak localization is described in terms of 
infinite impurity ladders called Cooperon- and Diffuson propagators [0, which satisfy the 
difi:usion equation. But the diffusive motion involves large changes in the direction of the 
particle due to many successive (possibly small angle) scatterings. Such a motion cannot be 
correctly described within the approximations inherent in higher-dimensional bosonization 
at the level of the Gaussian approximation |]5|,|6|, [TTl , p!2| (which amount (i) to the neglect of 
momentum-transfer between different patches, and (ii) to the local linearization of the energy 
dispersion 0), because in this case the electron trajectories are approximated by straight 
lines in the directions of the local velocities v". Obviously the weak localization effects must 
be contained in the corrections to the straight line approximation for the electron trajectory. 
Hence, the physics of weak localization will only emerge if we can generalize our approach 
such that it describes changes in the direction of the electron propagation due to successive 
scatterings. This can be achieved either by taking momentum-transfer between different 
patches (so-called around-the-corner processes) into account, or by considering patches with 
a finite curvature (i.e. by retaining the quadratic terms in the expansion of the energy 
dispersion at the Fermi surface) |]TB[. The equivalence of these procedures follows from the 
fact that any patch with a finite curvature can be subdivided into a larger number of ap- 
proximately fiat sub-patches, such that scattering from the original curved patch can also 
be described in terms of scattering from the set of coupled but fiat sub-patches. The latter 
point of view leads in the calculation of the Green's function for fixed background field 
to a system of coupled first order differential equations, which can always be transformed 
into a smaller system of differential equations with higher order derivatives. At present it 
is not clear whether it is possible to include the above corrections into our formalism in a 
simple approximate way, such that we do not loose the physics of weak localization. Obvious 
starting points might be to take around-the-corner processes between neighboring patches 
into account, or to retain second order spatial derivatives in the partial differential equation 
(0)- 

Let us now show that Eq.(|D reduces for static disorder to the usual perturbative result. 



In this case only the uOm = 0- component of the covariance function Cq is non-zero. For 

|qli 

simphcity let us assume that Cq has a separable form, Cq = S^^^oP'joe 9c , where |q|i = 
Y.i=i \Qi\ Then we obtain from Eq.(|^) for V ^ oo 

Qd..(r"v",r)~-^ , Kg.|»l , (11) 

where the inverse elastic mean free path is given 

by ^ = 2£ We conclude that at 
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large distances G°'(r, r) = Gg (r, r)e~ 2<^. In Fourier space this implies for |q| -C Qc 

G(k° + q,z^„) = -^ [ . , (12) 

where = is the inverse elastic lifetime. It is easy to see that this result agrees with the 
usual perturbative expression for the average self-energy in lowest order Born approximation, 
which is given by |l8l 



/mS(k) = ^ E e"^/mGo(k + q, -zO+) . (13) 

Because for kp only wave- vectors with |q| <^ kp contribute, it is allowed approximate 
the Green's function on the right-hand side of Eq.(|T3|) by its linearized form. Then we obtain 
— = 2/mE(k") = ( — ) in agreement with the above bosonization result. In Ref. 



|v-| 

T9| we have shown that a Debye- Waller factor which diverges at large distances stronger 
than logarithmic completely washes out any singularity in the momentum distribution n-^. 
Hence the average momentum distribution is for any finite disorder analytic at the Fermi 
surface. Of course, this is a well known result . It is also easy to understand why in one- 



dimensional interacting Fermi systems any finite disorder destroys the algebraic singularity 
of nk, which is one of the characteristics of a Luttinger liquid |2^. Recall that this algebraic 



singularity is due to the logarithmic divergence of Q^t(r, 0) for oo. At sufficiently large 
distances this weak logarithmic divergence is negligible compared with the linear divergence 
of g?,.(r,0). 

The case of a time-dependent random potential is more interesting. To calculate the 
average Green's function, we should specify the dynamic covariance function Cq. If we 



would like to model an underlying interacting many-body system in thermal equilibrium 
by a random system, then the form of Cg is completely determined by the nature of the 
interaction. In the case of the coupled electron-phonon system at high temperatures an ex- 
plicit microscopic calculation of Cg has been given by Girvin and Mahan |0. However, their 
identification of Cg with the parameters of the underlying many-body system is based on a 
perturbative calculation of the self-energy at high temperatures. Our functional bosoniza- 
tion approach allows us to relate the covariance function Cg of the random system in a much 
more direct and essentially non-perturbative way to the underlying many-body system. Ev- 
idently, the requirement that the average Green's function of the stochastic model should 
be identical with the Green's function of the many-body system without disorder is equiva- 
lent with Q2is{r,T) = Q°4(r,r), where Qrfi,(r,r) and (5°t(r,r) are given Eqs.(0) and (|10D. 



It immediately follows that the connection between the effective stochastic model and the 
interacting many-body system is given by the surprisingly simple relation 

C, = -fr^ . (14) 

For example, to describe longitudinal acoustic phonons with dispersion tUq that are coupled 
to the electrons via the Coulomb potential f^ = ^p-, we should choose the covariance 



function 21 



where 7 measures the strength of the electron-phonon coupling. We would like to emphasize 
that in spite of its apparent simplicity Eq.([T^) is a highly non-trivial result, because it is 
based on a controlled summation of the entire perturbation series of the many-body problem 
via bosonization. The crucial point is that bosonization produces an exponential resummation 
of the perturbation series, so that the effect of the interactions on the Green's function can be 
expressed exclusively in terms of a Debye- Waller factor Qf^f-{r,T) , which can then be directly 
compared with the Debye- Waller factor Q'^^g{r,T) due to disorder. 

Of course, the dynamic random potential could also be due to some non-equilibrium 
external forces, in which case the above identification with an underlying many-body system 
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is meaningless. As an example, let us consider a Gaussian white noise random potential, 

q|l 

with covariance given by Cg = C^e ic . Substituting this into Eq.(^ and taking the limit 
V, f3 oo, the integrations can be performed analytically, with the result 

iWr 

QLKv",r) = —-^ . . , , , , (16) 

where = 2^ (^) • Because the Debye- Waller factor vanishes at r = 0, we have 
G'°(r,0) = G'g(r, 0), so that the dynamic white noise random potential does not affect the 
momentum distribution nk. It is easy to see that this is an artifact of the white noise limit. 
Note, however, that for 00 

G"(r,r)~-5('^-i)(rl)— -— . (17) 

The above limit determines the density of states at the Fermi energy, which in turn can be 
expressed in terms of a renormalized effective mass. Hence, the random potential enhances 
the effective mass by a factor of e^/''^"'. This is intuitively clear: the electrons become 
heavier because they have to overcome the resistance of the random potential. The Fourier 
transformation of G°'{r, r) can be calculated exactly. For v" ■ q > the result is 

1 



G"(q,zcu) 



Wqc + zcj - v" ■ q 

(18) 



X < 1 H exp 

I — ■ q 



tLU — ■ q 

If we now analytically continue this expression to real frequencies by replacing iu u + iO'^, 
we encounter an essential singularity at = ■ q. We suspect that this is an artifact of 
the above simple choice of the random potential. 

In summary, we have shown that functional bosonization can be used to calculate the 
disorder averaged Green's function of electrons at high densities that are subject to a random 
potential with long-range spatial correlations. While in the static limit we have recovered 
the usual perturbative result, for time- dependent random potentials we have obtained a 
highly non-trivial expression for the averaged single-particle Green's function. One of our 
main results is Eq.(|l4D, which puts the description of an interacting many-body system via 



an effective stochastic model on a solid microscopic basis. Although our method describes 
disorder and interactions on equal footing, the corresponding contributions to the Debye- 
Waller factors are simply additive, so that interference terms containing weak localization 
effects do not appear. However, it might be possible to calculate the single-particle Green's 
function via higher-dimensional bosonization beyond the Gaussian approximation, taking 
momentum-transfer between different patches or the non-linearities in the energy dispersion 
approximately into account. For the density-density correlation function and the boson 
representation of the Hamiltonian such a calculation has been performed in Ref . |^ . In this 
case bosonization in > 1 could become a new powerful approach to disordered interacting 
Fermi systems in the diffusive regime, which can deal simultaneously with disorder and 
interactions and does not have the disadvantages of the replica trick. Note that the diffusive 
regime does not exist in = 1, so that the physics of weak localization can only be discussed 
within higher-dimensional bosonization. Work in this direction is in progress. 

The inclusion of transverse gauge fields is straightforward within our functional bosoniza- 
tion approach. The results of the present work can be combined with the method described 



in Ref. ||2^ to treat the problem of a random transverse gauge field with long-range spatial 
correlations in arbitrary dimensions. This problem is of experimental relevance in connection 
with the quantum Hall effect It should be mentioned, however, that the applicability 
of higher-dimensional bosonization to transverse gauge fields is controversial ||2^. We hope 
to address this issue in future publications. 

I would like to thank G. Castilla for carefully reading the manuscript and suggesting 
many useful modifications. I am also grateful to Kurt Schonhammer and Volker Meden for 
discussions. 
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